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Abstract 

One of severe limits of graphene nanoribbons (GNRs) in future applications is that zigzag GNRs 
(ZGNRs) are gapless, so cannot be used in field effect transistors (FETs). In this paper, using tight- 
binding approach and first principles method, we derived and proved a general edge (boundary) 
condition for the opening of a significant bandgap in ZGNRs with defective edge structures. The 
proposed semiconducting GNRs have some interesting properties including the one that they can 
be embedded and integrated in a large piece of graphene without the need of completely cutting 
them out. We also demonstrated a new type of high-performance all-ZGNR FET. 



1 



I. INTRODUCTION 



Graphene has attracted intensive research efforts due to its unique electronic and me- 
chanical properties.-"- A recent experiment demonstrated a beautiful technique in fabri- 
cating graphene nanoribbons (GNRs) with atomically precise edges, suggesting the great 
potential of GNRs in future applications of graphene-based high performance electronics.- 
Theoretical calculations showed that only two thirds of armchair graphene nanoribbons 
(AGNRs) with different widths are semiconducting, while zigzag graphene nanoribbons 
(ZGNRs) are gapless due to localized edge states at the Fermi level.'^ It was also theoret- 
ically demonstrated that these bandgap-closing edge states survive in zigzag GNRs with 
a mixture of zigzag and armchair sites at boundaries. Another recent theoretical work 
proved that the confinement by minimal boundaries generally does not produce an insu- 
lating GNR except for the armchair case.- This theoretically predicted edge or orientation 
and width dependence of bandgap opening in GNRs provide serious limits in real appli- 
cations of GNR-based electronic devices: First, gapless ZGNRs cannot be used in FETs, 
and second, the precise control of the width of AGNRs is required. 

In this paper, using the tight-binding approach and the first principles method based on 
density functional theory (DFT), we derived and proved that when the number of A-site 
defects equals to that of B-site defects at each boundary (A, B denotes two sublattices 
of graphene), localized edge states in GNRs will be eliminated, and then a bandgap that 
is inversely proportional to the ribbon width will generally be open. We then showed 
that ZGNRs with defective boundaries that satisfy the bandgap opening conditions can be 
embedded and integrated in a large piece of graphene, which may have implications for the 
future design of graphene-based integrated circuits. At last, we demonstrate a new type of 
field effect transistor completely made of ZGNRS. Note that in all previously theoretically 
proposed GNR-based transistors,- the AGNR is indispensable. It is worthy mentioning 
here that since the long-range magnetic order is not stable in one-dimensional systems 
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under finite temprature, we therefore stick to the non-magnetic case in this study. 

II. RESULTS AND DISCUSSION 

We first focus on a ZGNR with an edge structure with defects as shown in Fig.[TJa). The 
edge structure can be specified by a quadruple of segment lengths in unit of the graphene 
lattice constant (a = 2.46 A), (Nb^i, Naa, Nb,2, Na,2)- So the number of two-coordinated 
carbon atoms at edge belong to A (B) sublattice, Na (Nb), equals to Naa + Na,2 (Nb,i + 
Nb,2)- The electronic structure was calculated using the tight-binding approach. Only 
the nearest-neighbor hopping energy (-2.7 eV) was taken into account. It is well known 
that there exist edge states with £ = for a semi-infinite graphene with a zigzag edge 
if 27r/3 < kyU ^ n. The edge state entirely localizes at edge for kyU = n, and otherwise 
decays exponentially away from the edge.^ When two zigzag edges form a ZGNR and 
the edge states from both edges interact with each other, the edge states still degenerate at 
kya = n, while a small gap that decreases exponentially with the nanoribbon width opens 
elsewhere. The band structure of a perfect ZGNR in a supercell corresponding to (9, 0, 
0, 0) is reproduced in Fig. [lib). Due to band folding, there are six bands (marked as red) 
corresponding to 27r/3 < kyU ^ n. We find that these bands originating from edge states 
are removed and thus an energy gap opens if Na = Nb- An example of the band structure 
corresponding to a (3, 3, 3, 3) edge structure is shown in Fig. [He), and the squared wave 
functions in the inset clearly indicate they are extended states with the form, sm{knx), 
which has more nodes for larger energies. Therefore, the energy of the valence band 
maximum (VBM), hence the energy gap, is inversely proportional to the nanoribbon width. 
For the case of either Na < Nb or A^^ > A^s, some of the edge states will remain as shown 
in Figs.[IId) and[2e), though the perfectly localized edge state is destroyed. The squared 
wave functions in the insets of Figs.[2d) and[Ue) show the edge states are exponentially 
decaying away from the edge, so the energy gap due to the interaction between the states at 
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opposite edges also decreases exponentially with respect to the nanoribbon width. These 
two distinct behaviors of the bandgap variation as a function of the nanoribbon width can 
be seen in Fig.|2ta). The variation of bandgaps with respect to the characteristic length of 
the edge structure is shown in Fig. l2lb). The possibility to tune bandgaps with different 
edge structures on the same nanoribbon might provide useful implication in the design of 
nanoribbon-based electronic devices. 

The above mentioned condition for the elimination of edges states that leads to bandgap 
opening can be understood from the following arguments. Considering a semi-infinite 
graphene sheet with (A^|) two-coordinated carbon atoms at the edge and (N^) three- 
coordinated carbon atoms in the bulk belonging to A (B) sublattice, we have the following 
equation by the conservation of coordinate numbers, 

2Ni + 3N^ = 2N'B + 3N^. 

Since the total number of carbon atoms in each sublattice is = ^a(s) '^^A{By '•^^ 

above equation can be rewritten as 3(A'^ ~^b) ~ ^a~^b^ which means the relation be- 
tween A'^ and A^^ is the same as that between and N^. The band structure of graphene 
in the tight-binding approximation is calculated by 

£VA(r) =?[v|/fi(r)+\|/fi(r-Ri)+Vg(r-R2)] (1) 
E\\fB{r)=t[\\fA{r)+\\fA{r + Ri)+\\fA{r + R2)], (2) 

where t is the hopping energy, V|/A(r) and V|/g(r) are the wave functions on A and B atoms 
belonging to the same unit cell at a discrete coordinate r, and Ri and R2 are graphene 
lattice vectors as shown in Fig. [Ha). For E = 0, Eqs. [T]and|2]are decoupled. There are 
N'^ equations with A^^ unknowns for [T] and A^^ equations with A^^ unknowns for |2l So if 
N'^ > Ng, \T\ will have no solution while [2] will have solutions on A sublattice. Similar 
conclusion will arrive for N'^ < ^b- The fact that the wave function will reside on the 
sublattice with more atoms can be observed by comparison of the insets in Figs. [IJd) 
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and [He). If A^^ = Ng and the break of symmetry leads to no linear dependence among 
equations, then Eqs. [T] and [2] will have only zero solution, which is not admissible and 
results in the elimination of localized states with E = that decay exponentially from the 
edge. From the above arguments, it can be seen that this boundary condition is a necessary 
condition for opening bandgaps in GNRs that scale inversely with the width instead of 
exponentially. Our test calculations also show that this boundary condition of bandgap 
opening also applies for many other nanoribbons with general orientations and defective 
edge structures. 

In practice, disorders are inevitable in the fabrication of GNRs that involves in the 
lithographic patterning and etching. We therefore considered a disordered edge structure 
as shown in Fig. IH A ZGNR with one (3, 3, 2, 3) unit (as shaded in the figure) plus 10 
(3, 3, 3, 3) units in one supercell. Here, the (3, 3, 2, 3) unit can be treated as a disorder in 
the (3, 3, 3, 3) edge structure. Our tight-binding calculations show that the disorder will 
induce localized states inside the bandgap (as shown in the figure) that may not contribute 
to transport, and compared to the edge structure without the disorder, the energy gap be- 
tween extended states in the disordered structure is bigger than the bandgap of the perfect 
one. In real experiments, statistically speaking, the gap-opening boundary condition is 
always satisfied, so our findings presented here may have implications for recent experi- 
ments showing that the transport gap of GNRs is inversely proportional to the width, and 
independent on the orientations or edges of GNRs.— 

An interesting property of nanoribbons with a bandgap-opening edge structure is that 
if a wide nanoribbon is joined with a narrow nanoribbon, the electronic structure of the 
wide nanoribbon near the Fermi energy is not altered with electrons still confined in the 
wide nanoribbon. An example of a nanoribbon with a width of L = 12\/3a and a (3, 3, 
3, 3) edge structure joined with a nanoribbon with the same edge structure and a width 
of L = 5v^a is shown in Fig. [3ta). The band structure of the compound system near the 
Fermi energy (the conduction and valence band) in Fig. |3lb) is almost the same as that 
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of the stand-alone nanoribbon shown in Fig. [He). The charge distribution of the state at 
VBM in Fig.[3ta) and the local density of states in Fig.[3tc) clearly indicate that the wave 
function is only localized in the wide nanoribbon. The confinement can be understood 
from the bandgap difference of two nanoribbons with different widths. Note that the in- 
tegrated GNRs discussed here can be also regarded as a special type of graphene antidot 
lattice structures proposed earlier.— This property makes it possible to fabricate indi- 
vidual nanoribbon-based electronic devices by patterning rows of holes in a large piece 
of graphene avoiding complete cutout and glued together, which might be beneficial for 
the integration of future graphene-based electric circuits. On experimental side, the pat- 
terned graphene nanostructures discussed here can be obtained experimentally by using 
techniques such as templated self-assembly of block copolymers'^ or direct writing using 
a helium ion beam.— The periodically patterned structure may be formed by first form- 
ing resist patterns on the graphene sheet followed by templated self-assembly of block 
copolymers in the region where the resist have been removed and etching of graphene by 
using the copolymer patterns as the mask. On the other hand, the random patterns can be 
formed by direct writing using a helium ion beam. Prior to the lithography processes, an 
alignment mark may be formed on the wafer by using an appropriate graphene edge as a 
reference so as to align the patterns in specific directions with respect to the underlying 
graphene lattice structure. 

At last, we show a FET completely made of ZGNRs as shown in Fig. [51 In all previous 
theoretically proposed GNR-based FETs, the AGNR is an indispensable component due 
to the fact that pure ZGNR is metallic. Here, the proposed FET consists of two pure 
ZGNR electrodes (left and right), and a ZGNR with a defective (2, 2, 2, 2) edge structure. 
The transport calculations were done using a first principles approach combining the non- 
equilibrium Green's function's techniques and DFT via the ATK code.—"— In the inset 
of the figure, the current-voltage (I-V) curve is shown for the zero gate voltage. The bias 
range of the zero current comes from the bandgap of the defective ZGNR in the center, 
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confirming the bandgap opening condition we derived from the tight binding approach. 
The currents as a function of gate voltage for different bias voltages suggest that the on-off 
ratio of this proposed FET is bigger than 1000. Compared to the previously proposed all- 
GNR based FET that used two ZGNRs and one AGNR,^ the FET suggested here has two 
obvious advantages: First, the complicated contacts between differently orientated AGNR 
and ZGNRs are avoided. Second, the precise control of ribbon width is not required. 

III. CONCLUSION 

In conclusion, using the tight-binding approach, we derived a general boundary condi- 
tion for the band gap opening in the ZGNRs with defective edges: When the number of 
A-site defects equals to that of B-site defects, the ZGNRs are semiconducting. We further 
showed that the semiconducting ZGNRs generated this way can be integrated in a large 
piece of graphene by correctly patterning holes, which may be useful for the future large- 
scale integration of GNR-based devices. At last, we demonstrated using first principles 
calculations a high-performance FET completely made of ZGNRs. Results presented in 
this paper may be used to explain the recent experimental measurements showing that the 
transport gap always exists independent of the crystallographic orientations of GNRs. We 
expect these findings to provide impetus for new experiments as well motivations for new 
ideas in designing ZGNR-based electronic devices. 
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FIG. 1 : (a) The lattice structure of a zigzag graphene nanoribbon with periodic edge structures. 
The unit cell is indicated by the dashed Une. The edge structure is denoted by a quadruple of 
(Nb,i, Na,i, Nb,2, Na,2), each number of which corresponds to the segment length in unit of the 
graphene lattice constant, a. Other parameters of the system are the nanoribbon width (L) and the 
translational vector (T). Carbon atoms belonging to different sublattices at edge are designated red 
(A) and blue (B) colors, (b-e) The band structures of nanoribbons with different edge structures 
and the same width L = 12\/3a. In the inset of (c) are shown the squared wave functions along the 
dash-dotted Une in (a) for different states at ky = 0. The squared wave functions corresponding to 
the valence band maximums are also plotted in the insets of (d) and (e). The radii of filled discs are 
proportional to /?(logjQ |\|/(r)p +4), where R{x) is the ramp function, and the color is determined 
by the sign of real part of \|/(r). 
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FIG. 2: (a) The variation of bandgaps as a function of the nanoribbon width (L) for different edge 
structures, (b) The variation of bandgaps as a function of the segment length (Na) of a (N, N, N, 
N) edge structure for different widths. 
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FIG. 3: (a) The unit cell (dashed line) of a wide nanoribbon (W) joined with a narrow nanoribbon 
(A^) having the same edge structure. The nanoribbon W between two dash-dotted lines is identical 
to the nanoribbon shown in Fig. [IJa). (b) The band structure of the system in (a). The squared 
wave function plotted in (a) corresponds to the state indicated by an arrow, (c) The corresponding 
local density of states in W and N as shown in (a). 
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FIG. 4: The comparison between band structures of graphene nanoribbons with a pure (3, 3, 3, 3) 
edge structure (dashed hne) and with one (3, 3, 2, 3) unit (as shaded) plus ten (3, 3, 3, 3) units in a 
supercell (solid line). The squared wave function (only part of the supercell is depicted) shown in 
the left panel suggests that the corresponding state (as arrowed) is localized around the defective 
(3, 3, 2, 3) unit. 
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FIG. 5 : (a) The configuration for the transport calculation, (b) The variation of current as a function 
of gate voltage for different bias voltages. The inset shows the bias voltage dependence of current 
for zero gate voltage. 
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